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Abstract 

Background: The molecular bases of mammalian pancreatic a cells higher resistance than (3 to proinflammatory 
cytokines are very poorly defined. MicroRNAs are master regulators of cell networks, but only scanty data are 
available on their transcriptome in these cells and its alterations in diabetes mellitus. 

Results: Through high-throughput real-time PGR, we analyzed the steady state microRNA transcriptome of murine 
pancreatic a (aTGl-6) and (3 (PTGl) cells: their comparison demonstrated significant differences. We also 
characterized the alterations of aTGl-6 cells microRNA transcriptome after treatment with proinflammatory 
cytokines. We focused our study on two microRNAs, miR-296-3p and miR-298-5p, which were: (1) specifically 
expressed at steady state in aTGl-6, but not in PTGl or INS-1 cells; (2) significantly downregulated in aTGl-6 cells 
after treatment with cytokines in comparison to untreated controls. These microRNAs share more targets than 
expected by chance and were co-expressed in aTGl-6 during a 6-48 h time course treatment with cytokines. The 
genes encoding them are physically clustered in the murine and human genome. By exploiting specific microRNA 
mimics, we demonstrated that experimental upregulation of miR-296-3p and miR-298-5p raised the propensity to 
apoptosis of transfected and cytokine-treated aTGl-6 cells with respect to aTGl-6 cells, treated with cytokines after 
transfection with scramble molecules. Both microRNAs control the expression of IGF1R|3, its downstream targets 
phospho-IRS-1 and phospho-ERK, and TNFa. Our computational analysis suggests that MAFB (a transcription factor 
exclusively expressed in pancreatic a cells within adult rodent islets of Langerhans) controls the expression of 
miR-296-3p and miR-298-5p. 
(Continued on next page) 
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(Continued from previous page) 

Conclusions: Altogether, high-throughput microRNA profiling, functional analysis with synthetic mimics and 
molecular characterization of modulated pathways strongly suggest that specific downregulation of miR-296-3p 
and miR-298-5p, coupled to upregulation of their targets as IGF1R|3 and TNFa, is a major determinant of 
mammalian pancreatic a cells resistance to apoptosis induction by cytokines. 

Keywords: Mammalian pancreatic a and (3 cells, microRNA transcriptome. Proinflammatory cytokines, Apoptosis, 
Cellular networks. Diabetes mellitus 



Background 

Insulitis is an inflammation of pancreatic Langerhans islets, 
which is known to precede the onset of diabetes mellitus 
(DM) [1,2]. Post-insulitis decrease of pancreatic |3 cells is a 
major hallmark of both type 1 (TIDM) and type 2 diabetes 
mellitus (T2DM) [3-5]. It is associated to a cells dysfunc- 
tions and high glucagon secretion, which contribute to 
chronic hyperglycemia and ensuing clinical outcomes in 
DM [6,7]. Accordingly, efforts to increase our knowledge 
on the biomolecular mechanisms regulating proliferation, 
physiopathological functions, and apoptosis of a cells will 
likely result in improved medical approaches to the disease. 
The pathways leading to mammalian p cell apoptosis, 
induced by proinflammatory cytokines, are complex but 
sufficiently known [8,9]. On the contrary, the molecular 
bases of a cells higher resistance than |3 to the same cues 
are still largely uncharacterized [10-12]. Pancreatic a and |3 
cells share common endocrine precursors (Ngn3^ cells), 
but upon differentiation they respond differently to external 
stimuli as proinflammatory cytokines, glucagon-like peptide 
1 (GLPl), and inhibitors of dipeptidyl peptidase-4 (DPP-4) 
[12-14]. MicroRNAs (miRNAs) are small (19-25 nucleo- 
tides) noncoding RNAs, which control gene expression 
mainly at the post-transcriptional level: miRNAs have been 
shown to be master regulators of cell networks [15,16]. 
Many miRNAs are involved in DM physiopathology: for in- 
stance, miR-375 controls |3 cell mass and insulin secretion 
[17]; miR-21 and miR-146a are involved in |3-cell apoptosis 
[18]; miR-30d promotes insulin synthesis and protects |3 
cells from damage by proinflammatory cytokines [19]; miR- 
126 is responsible for impaired angiogenetic signaling in 
DM patients [20]. However, there are no high- throughput 
(HT) published data on a cells miRNA transcriptome at 
steady state, its comparison with |3 cells steady state miRNA 
transcriptome, and its alterations following treatment with 
cytokines. To discover new physiopathologic regulatory 
mechanisms, specific of cytokine-treated pancreatic cells, 
we profiled global miRNAs expression in two mouse pan- 
creatic cell lines (aTCl-6 and |3TC1) at steady state and in 
aTCl-6 after treatment with proinflammatory cytokines 
IFN-y, IL-1|3, TNF-a. This analysis brought to the identifl- 
cation of: (1) differentially expressed (DE) miRNAs between 
aTCl-6 and pTCl at steady state; (2) DE miRNAs between 



aTCl-6 treated with cytokines and matched untreated con- 
trols. To insert our data in the appropriate biological con- 
text, we reconstructed the molecular networks regulated at 
steady state by DE miRNAs in a and |3 cells. Our results 
demonstrate that these networks are remarkably different: 
very likely, they increase a cells resistance to apoptosis and 
negatively interfere with pancreatic |3 cells viability after ex- 
posure to cytokines. 

Results 

Apoptosis induction by cytokines in aTC1-6 and piCI 
cells 

To verify the differential response of aTCl-6 and pTCl 
cells to treatment with cytokines, we assayed their propen- 
sity to undergo apoptosis during a time course experiment 
(6-24-48 h after treatment, PT). The number of apoptotic 
aTCl-6 cells did not significantly vary with respect to con- 
trols for the entire time course (Figure lA), demonstrating 
that in our system a cells appear to withstand apoptosis 
induction by cytokines; on the contrary, |3TC1 cells are 
clearly more susceptible to treatment (Figure IB). 

Steady state miRNA transcriptome profiles of aTC1-6 and 
PTCl cells 

To identify miRNAs potentially responsible for aTCl-6 
resistance to cytokines, we analyzed the steady state 
miRNA transcriptome of aTCl-6 and pTCl ceUs. At 
steady state, 23 miRNAs were exclusively expressed in 
aTCl-6 (a-miRNAs), while 26 were expressed only in 
pTCl (|3-miRNAs); 50 miRNAs resulted significantly 
more expressed in aTCl-6 than in |3TC1, whereas 74 
were significantly more expressed in pTCl compared to 
aTCl-6 (Limma test, Benjamini-Hochberg adjusted 
p-values < 0.05) (see Additional file 1). Their assignment 
to specific famflies is shown in Additional files 2 and 3. 

Treatment with cytokines alters aTC1-6 miRNA 
transcriptome profiles 

After treatment with cytokines for 48 h, 3 miRNAs 
(miR-146a, miR-203, miR-298-5p) were significantly dif- 
ferentially expressed in aTCl-6 as compared to matched 
untreated controls (Limma test, Benjamini-Hochberg 
adjusted p-values < 0.05) (see Additional file 4). Their 
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Figure 1 Apoptosis of aTC1-6 and PTC1 after treatment with cytol<ines. (A) Annexin V flow cytometric analysis of apoptosis in aTCl-6 
treated with IFN-y, IL-1(3, TNF-a for 6, 24, 48 h and in matched untreated controls. Values represent the percentage of apoptotic cells. Data are 
presented as mean ± S.D. of three independent experiments (n = 3). (B) Microphotographs of (3TC1 at steady state (left) and 24 h PT (right). Cell 
shrinkage and irregular morphology are evident in (3TC1 after 24 h of treatment. Representative pictures are shown from three independent 
experiments (n = 3). 



assignment to specific families is shown in Additional 
file 5. 

MiRNAs 296-3p and 298-5p are reliable candidates for 
involvement in aTC1-6 higher resistance than (3 cells to 
apoptosis induction by cytokines 

To identif)^ miRNAs whose functions could explain the 
differential response to cytokines of pancreatic a and p 
cells, we specifically focused our attention on miR-296- 
3p and miR-298-5p. Both miRNAs are expressed at 
steady state only in aTCl-6, but are not synthesized in 
pTCl (Figure 2A); both are part of the imprinted Gnas/ 
GNAS clusters in mice and humans and share more tar- 
gets than expected by chance (44 versus 33, respectively; 
p = 0.0485, /^-test), even though their seed regions are 
different: this suggests common regulatory functions. 
Single TaqMan gene expression assays (STAs) during a 
time course analysis (6-12-24-48 h) showed a highly sig- 
nificant downregulation of miR-296-3p in aTCl-6 cells 
at 24 and 48 h PT, compared to matched untreated con- 
trols (Student s ^-test, Bonferroni adjusted p-value < 
0.01). Mir-298-5p resulted significantly downregulated 
starting at 12 h PT and reached a highly significant 
downregulation at 48 h PT, compared to matched un- 
treated controls: Students ^-test, Bonferroni adjusted p- 
values, were <0.05 and <0.01, respectively (Figure 2B). 
Both miRNAs were co-expressed in aTCl-6 cells through- 
out the entire experimental time course (r-value = 0.88, p- 
value = 1.15e-08, Pearsons correlation test) (see Additional 
file 6). STAs confirmed that they are not expressed either 
in |3TC1 or INS-1 cells. 



Genomics of genes encoding miR-296-3p, miR-298-5p, 
Nespas and identification of upstream CpG islands 

Genes encoding miRNAs 296-3p and 298-5p are clustered 
in a genomic region, which also comprises the gene for 
the noncoding transcript Nespas and is imprinted in mice 
and humans [21]. Sequences of mature miR-296-3p are 
100% conserved between rodents and humans, whereas 
those of miR-298-5p are 74% identical. Analysis of this re- 
gion through UCSC browser revealed the presence of two 
clusters of CpG islands: (i) one comprises two CpG 
islands, from 17.5 to 18.8 kb upstream the first nucleotide 
of pre-miR-296, and is located 9.2 and 10 Kb downstream 
Nespas transcription start site (TSS); (ii) the other is made 
of three CpG islands from 30.1 to 33.6 Kb upstream the 
first nucleotide of pre-miR-296 and is located 1.8-5 kb up- 
stream Nespas TSS (see Additional file 7). Matlnspector 
revealed a putative Nespas promoter located 500 bp 
upstream- 100 bp downstream its TSS. 

aTC1-6 transfection with mimics of miR-296-3p and miR- 
298-5p increases apoptosis levels induced by cytokines 

To precisely define miR-296-3p and miR-298-5p biological 
functions, we compared apoptosis levels of aTCl-6 cells, 
transfected with either one or both miRNA mimics and 
treated with cytokines for 6, 24, 48 h after transfection 
(AT), with those of scramble-transfected aTCl-6 cells 
treated with cytokines by following a similar protocol. The 
percentage of apoptotic aTCl-6 cells after transfection 
with mimics of miR-296-3p was comparable to scramble- 
transfected controls during the entire time course treat- 
ment. On the other hand, transfection with mimics of 
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Figure 2 Expression of miR-296-3p and miR-298-5p in aTC1-6 and piCI. (A) Real-time PGR amplification plot of miR-296-3p and miR-298-5p 
in aTCl-6 and (3TC1 at steady state. Note that both miRNAs are not expressed in (3TC1. Plot image is representative of three independent 
experiments (n = 3). (B) Gene expression fold changes of miR-296-3p and miR-298-5p in aTCl-6 at 6, 12, 24, 48 h PT, relative to matched 
untreated cells. MiRNAs expression was measured by quantitative real-time PGR. MiR-26a was used as endogenous control. Data are reported as 
mean ± S.D. of three independent experiments (n = 3). * p-value < 0.05; '''' p-value < 0.01 (Student's t-test, Bonferroni correction). 




miR-298-5p or of both miR-296-3p and miR-298-5p 
increased in a highly significant manner the number of 
aTCl-6 apoptotic cells, compared to matched scramble- 
transfected controls: 1.5 and more than 2.5 folds at 24 h 
PT, respectively (Tukey HSD post-hoc one-way ANOVA 
test, p-value <0.01) (Figure 3). 

Identification of miR-296-3p and miR-298-5p targets 

To characterize the networks regulated by miRNAs 296- 
3p and 298-5p, we computationally searched their vali- 
dated and predicted targets. We identified 1 validated 



target of miR-296-3p; 5 validated targets of miR-298-5p; 
207 predicted targets of miR-296-3p; 707 predicted tar- 
gets of miR-298-5p. We focused our attention on 7 tar- 
gets of miR-296-3p, 4 of miR-298-5p, 2 common to both 
miRNAs: they were chosen according to their involve- 
ment in apoptosis, cell cycle progression, cell differenti- 
ation and hormone secretion (see Additional file 8). 
Among them, Bcl2, Ccna2, Irs2, Nr4a2 are transcription- 
ally regulated by CREBl, which is a validated target of 
miR-296-3p [22]; Tnf^xvd Vdr are validated targets of miR- 
298-5p [23,24]. All other targets were computationally 
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Figure 3 Upregulation of miR-296-3p and miR-298-5p reduces aTC1-6 resistance to apoptosis induced by cytolcines. Annexin V flow 

cytometric analysis of apoptosis in aTCl-6 transiently transfected with scrambled oligonucleotides (NC, Negative Control), mimics of miR-296-3p, 

miR-298-5p, a mix of both, and treated with IFN-y, IL-1 (3, TNF-a for 6, 24, 48 h. The y-axis represents the percentage of apoptotic cells; the x-axis 

represents the four experimental conditions [(i)scramble-transfected cells; (ii) cells transfected with mimics of miR-296-3p; (iii) cells transfected 

with mimics of miR-298-5p; (iv) cells transfected with a mix of both] assayed at the three time-points. Data are presented as mean + S.D. of three 

independent experiments (n = 3). All the possible pairwise comparisons were performed among the four different experimental conditions within 

each time point: significant differences have been assessed through Tukey HSD post-hoc one-way ANOVA test p-value < 0.01). At 24 h PT, 

aTCl-6 transfected with mimics of miR-298-5p show a highly significant increase of the number of apoptotic cells with respect to scramble- 

transfected control; at the same time point, in aTCl-6 transfected with mimics of both miR-296-3p and miR-298-5p a highly significant increase of 

the number of apoptotic cells is detected with respect to all the other experimental conditions. 
V J 



predicted, including IGFlRp that we validated through 
western analysis (see later). 

Modulation of miR-296-3p and miR-298-5p also alters 
expression of their targets 

To verify whether in vitro modulation of miR-296-3p 
and miR-298-5p affected the expression of their targets, 



we performed transient transfection experiments of 
aTCl-6 cells with their mimics. Transfection efficiency 
at 24 h AT was in all cases higher than 90%. In aTCl 
transfected with mimics of miR-296-3p or miR-298-5p, 
real-time PGR showed altered expression of different 
genes with respect to scramble-transfected cells, includ- 
ing Igflr, Tnf, Vdr (Figure 4). 



1.6 



O 

^ 0.5 



§ 0 



-0.5 



-2.1 



_[L_ 




■ aTC1-6 24HCyt ■ aTC1-6 296-3p 24H No Cyt ■ aTC1-6 298-5p 24H No Cyt 
□ aTC1-6 48HCyt □ aTC1-6 296-3p 48H No Cyt □ aTC1-6 298-5p 48H No Cyt 



T 



Figure 4 Modulation of miR-296-3p and miR-298-5p alters expression of their targets. Bar graph showing changes in gene expression of a 
selected set of miR-296-3p and miR-298-5p targets for each of three different experimental conditions: (i) aTCl-6 treated with cytokines with 
respect to matched untreated control cells at the time points 24 and 48 h; (ii) untreated aTCl-6 transfected with mimics of miR-296-3p with 
respect to scramble-transfected control cells at the time points 24 and 48 h; (iii) untreated aTCl-6 transfected with mimics of miR-298-5p with 



respect to scramble-transfected control cells at the time points 24 and 48 h. Data are reported as LOG of 2^ 
as endogenous controls to normalize real-time PCR data. 



values. Hprt and Ppia were used 
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Expression of IGF1 RP and TNFa is controlled by miR-296- 
3p and miR-298-5p in aTC1-6 cells 

In aTCl-6 treated with cytokines for 24 h, levels of pro- 
teins IGFlRp and TNFa increased 1.5 and 17 folds with 
respect to untreated controls, respectively (Figures 5A 
and 5B, right panels), while miRNAs 296-3p and 298-5p 
decreased of 3.7 and 2.2 folds with respect to untreated 
controls, respectively (Figure 2B). We further demon- 
strated through western analysis that protein IGFlRp 
decreased about 1.4 folds in steady-state aTCl-6 after 
transfection with mimics of miR-296-3p or miR-298-5p, 
as compared to scramble-transfected controls; this de- 
crease was higher than 3 folds when cells were trans - 
fected with both mimics (Figure 5 A, left panel). The 
decrease of IGFlRp was not detectable in aTCl-6 trans- 
fected with either mimic of miR-296-3p or miR-298-5p 
and then treated with cytokines for 24 h. In aTCl-6 
treated with cytokines for 24 h, after transfection with 
mimics of both miRNAs 296-3p and 298-5p, it was in- 
stead lower than at steady state (1.4 folds). This could 
also be due to the decrease of miRNAs 296-3p and 298- 
5p following treatment with cytokines (Figures 2B and 
5A, middle panel). Following transfection of aTCl-6 



with either mimics of miR-296-3p or miR-298-5p, TNFa 
protein was about 1.2 folds less expressed with respect to 
scramble-transfected controls; this decrease was more 
pronounced (1.6 folds) when cells were transfected with 
both mimics (Figure 5B, left panel). Each of the two 
miRNA mimics decreased the expression of TNFa more 
than 2 folds after 24 h of treatment with cytokines, but the 
silencing effect was not potentiated when the two mimics 
were transfected together (Figure 5B, middle panel). 

Activation of IRS-1 and ERK-1 is also under control of 
miR-296-3p and miR-298-5p 

To verify if the IGFIR signaling pathway was controlled 
by miR-296-3p and miR-298-5p through IGFlRp, we 
assayed the phosphorylation levels of IRS-1 and ERK (two 
markers downstream the IGF-1 receptor) in aTCl-6 
transfected with mimics of each one or both miRNAs; this 
analysis was performed both after treatment with cyto- 
kines for 24 h or on untreated cells, exploiting as controls 
matched scramble-transfected aTCl-6 cells. The expres- 
sion of phospho-IRS-1 didn t change in untreated aTCl-6 
transfected with mimics of each miRNA alone, while it 
decreased about 1.5 folds in aTCl-6 transfected with 
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Figure 5 Expression of IGF1RP and TNFa proteins is regulated by miR-296-3p and miR-298-5p in aTC1-6. (A) Western analysis of IGF1R(3 
in (1) untreated aTCl-6 transfected for 24 h with (i) scramble molecules (NC); (ii) mimics of miR-296-3p; (iii) mimics of miR-298-5p; (iv) mimics of 
both miR-296-3p and miR-298-5p (left); (2) aTCl-6 transfected for 24 h with (i) scramble molecules (NC); (ii) mimics of miR-296-3p; (iii) mimics of 
miR-298-5p; (iv) mimics of both miR-296-3p and miR-298-5p and treated with cytokines for further 24 h (middle); (3) aTCl-6 treated with 
cytokines for 24 h and their matched untreated controls (right). (B) Western analysis of TNF-a performed in the same experimental conditions as 
(A). |3-Actin signal was used to normalize the data. Numbers below Actin blots represent fold change expression values relative to 
matched controls. 
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mimics of both miR-296-3p and miR-298-5p (Figure 6A, 
left panel). By using the same controls, we detected a 
slight decrease of phospho-IRS-1 (1.2 folds) in aTCl-6 
transfected with mimics of both miR-296-3p and miR- 
298-5p and treated with cytokines (Figure 6A, right panel). 
Interestingly, also the activation of ERK-1 appears to be 
regulated by miR-296-3p and miR-298-5p: in the absence 
of treatment with cytokines, aTCl-6 cells transfected with 
mimics of miR-296-3p showed levels of phospho-ERK-1 
(Thr202) similar to those found in scramble-transfected 
aTCl-6 cells; the transfection with mimics of miR-298-5p 
or of both miRNAs led instead to a decrease of the protein 
(1.2 and 1.3 folds, respectively) (Figure 6B, left panel). 
Under cytokine treatment, the amount of phospho-ERK-1 
(Thr202) did not change in aTCl-6 transfected with either 
mimics of miR-296-3p or of miR-298-5p alone, compared 
to scramble-transfected cells; in aTCl-6 simultaneously 
transfected with mimics of both miRNAs, we detected a de- 
crease of 1.7 folds of phospho-ERK-1 (Thr202) levels, com- 
pared to scramble-transfected cells (Figure 6B, right panel). 

Identification of TFs regulating DE miRNAs 

In the genomic region comprising the genes encoding 
miR-296-3p, miR-298-5p and Nespas, Matlnspector iden- 
tified Transcription Factor Binding Sites (TFBS) for sixty 
seven Transcription Factors (TFs). Four of them (HMX2, 
HNF4A, LEFl, MAFB) are known to be expressed in the 
islets of Langerhans, and MAFB is known to be expressed 
only in rodent islet a cells within adult pancreas [25]; the 
presence of binding sites for this TF within the promoter 
of the genes encoding miR-296-3p and miR-298-5p sug- 
gests that it may regulate the expression of both miRNAs. 
Multi Experiment Matrix (MEM) showed a statistically 
highly significant negative correlation (p- values < 0.0001, 
Pearson Correlation test) between Maft) and Igflr. Interest- 
ingly, Microrna.org predicts two binding sites for mmu-miR- 
296-3p and mmu-miR-298-5p on Maft> mRNA 3' UTR. 

Differences between aTC1-6 and piCI steady state 
regulatory networks 

14 of 73 miRNAs, more abundantly or exclusively ex- 
pressed at steady state in aTCl-6 with respect to pTCl, 
have validated targets known to be expressed in mam- 
malian pancreatic a cells (see link to TlDbase in Materi- 
als and Methods). By using these data, we generated a 
network comprised of 560 nodes and 8222 edges (see 
Additional file 9): from this network, we inferred a sub- 
network made of 4 a-miRNAs for which validated tar- 
gets have been identified, consisting of 117 nodes and 
530 edges (see Additional file 10). In the case of pTCl 
cells, we found that 36 of 100 miRNAs more abundantly 
or exclusively expressed with respect to aTCl-6 have 
validated targets known to be expressed in pancreatic p 
cells (see link to TlDbase): a network of 439 nodes and 



2079 edges was generated from them (see Additional file 
11). The subnetwork made of 15 p-miRNAs consisted of 
107 nodes and 132 edges (see Additional file 12). The en- 
richment in specific biological processes was calculated in 
pancreatic aTCl-6 versus pTCl cells subnetworks (net- 
works made of genes regulated by miRNAs specifically 
expressed at steady state in aTCl-6 and pTCl cells, re- 
spectively) (Table 1). Interestingly, in the network of genes 
regulated by a-miRNAs, pancreatic a cell fate commit- 
ment is a biological process overrepresented with respect 
to the network regulated by p-miRNAs (Hypergeometric 
test, Benjamini-Hochberg adjusted p-value < 0.05). Within 
the network of genes regulated by a-miRNAs, insulin-like 
growth factor receptor signaling pathway is a biological 
process significantly enriched among the genes interacting 
with targets of miR-296-3p and miR-298-5p, with respect 
to the genes linked to the targets of the other miRNAs (p- 
value = 0.039, Fisher s exact test). 

Discussion 

The data reported in this paper and those from the litera- 
ture [12] confirm that mammalian pancreatic a cells are 
more resistant than p cells to apoptosis induced by cyto- 
kines (Figures lA and IB). Our results suggest that miR- 
296-3p and miR-298-5p play a pivotal role in determining 
this trait. Due to DM epidemic spreading [26], its molecular 
characterization has become an important biological and 
translational problem in BioMedicine [7]. Our HT and net- 
work analyses of miRNA transcriptome in a and p cells 
highlighted very different profiles at steady state (see Table 1, 
Additional file 1). Computational analysis of these data 
allowed us to identify a set of miRNAs, specifically synthe- 
sized in pancreatic a cells, which appear to also regulate 
pancreatic a cell commitment. Specific enrichment of path- 
ways, regulated by a-miRNAs within this biological process 
(Table 1), allowed the identification of miRNAs that could 
be involved in the acquisition of the trait resistance to 
cytokines-induced cell death upon pancreatic a cells differ- 
entiation. Among them, miR-296-3p and miR-298-5p stood 
out clearly as potentially critical nodes, responsible for a 
cells resistance to cytokine-induced cell death. In fact, be- 
sides being specifically expressed at steady state in a cells 
(aTCl-6), but not in p cells (pTCl or INS-1), both miRNAs 
were significantly co-expressed and downregulated in 
aTCl-6 during a time-course treatment with cytokines, 
compared to matched untreated controls (see Figures 2A 
and 2B, Additional file 6). Our subsequent functional assays 
demonstrated that aTCl-6 transfected with mimics of 
miR-296-3p and miR-298-5p became susceptible to apop- 
tosis when treated with the same cocktail of cytokines as 
compared to aTCl-6 transfected with scramble molecules. 
Although its biological effect was potentiated when also 
miR-296-3p was expressed, our results suggest that the role 
of miR-298-5p is more important in this process than that 
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Figure 6 Activation of IRS-1 and ERK-1 is under control of miR-296-3p and miR-298-5p in aTC1-6. (A) Western analysis of phospho-IRS-1 (Jyr612) 
in (1) untreated aTCl-6 transfected for 24 li witli (i) scramble molecules (NC); (ii) mimics of miR-296-3p; (iii) mimics of miR-298-5p; (iv) mimics of both 
miR-296-3p and miR-298-5p (left); (2) aTCl-6 transfected for 24 h with (i) scramble molecules (NC); (ii) mimics of miR-296-3p; (iii) mimics of miR-298-5p; (iv) 
mimics of both miR-296-3p and miR-298-5p and treated with cytokines for further 24 h (right). (B) Western analysis of phospho-ERK-1/2 ^hr202^yr204) 
performed in the same experimental conditions as in (A). Quantification of immunoblot signals was made by equalizing phospho-specific IRS-1 or Erkl/2 
band intensities to total IRS-1 or Erkl/2, respectively. The decrease in phosphorylation was normalized to the basal level of the control and reported in 
arbitrary units as fold decrease over basal value. Numbers below Actin blots represent fold change expression values relative to matched controls. 



of miR-296-3p (Figure 3): this would highlight the role of 
one or a few specific miR-298-5p targets. Through western 
analysis, we confirmed our computational prediction that 
IGFlRp and TNFa are common targets to both miRNAs 
and that miR-296-3p and miR-298-5p also control IRS-1 



Table 1 Top 5 biological processes significantly enriched 
in a-miRNA networks with respect to P-miRNA networks 



Biological process 


Adj-p-value 




(Benjamini-Hochberg) 


DNA replication 


3.15E-23 


Cell division 


2.51E-09 


Cellular response to metal ion 


4.61 E-05 


Positive regulation of mesenchymal cell 


5.24E-05 


proliferation 




Pancreatic a cell fate commitment 


2.23E-02 



Biological processes are ranked based on increasing adjusted p-values (from 
the most to the least significant biological process). 



and ERK-1 within the IGFIR signaling pathway. IGFIR is 
known to promote resistance to apoptosis by different 
mechanisms: it increases the levels of antiapoptotic pro- 
teins, as BCL2 and BCL-XL; it inactivates proapoptotic pro- 
teins, as BAD and CASP9; it stimulates mitogenic IRS-1/ 
MAPKs pathway [27-29]. Phospho-IRSl (Tyr612) and 
Phospho-ERK-1 (Thr202), which we demonstrated to be 
controlled by miR-296-3p and miR-298-5p, are known to 
regulate the response to insulin and to be involved in sur- 
vival and proliferation processes [30]. Decreased expression 
of mir-296-3p and miR-298-5p and the corresponding acti- 
vation of survival and proliferation signals, mediated by 
IGFIR and its downstream nodes (e.g., IRS-1 and ERK-1), 
may thus explain why aTCl-6 cells are resistant to death 
induction by cytokines (see Additional file 13). TNFa is 
known to be physiologically expressed by pancreatic endo- 
crine cells and to contribute to maintain islet homeostasis 
[5,31]. This protein is also known to decrease apoptosis levels 
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of pancreatic acinar cells during acute pancreatitis by stimu- 
lating the synthesis of antiapoptotic proteins [32]. 

Conclusions 

HT miRNA profiling data, functional analysis with syn- 
thetic mimics and molecular characterisation of modulated 
pathways strongly suggest that specific downregulation of 
miR-296-3p and miR-298-5p in pancreatic a cells, coupled 
to upregulation of their targets as IGFlRp and TNFa and 
activation of the corresponding signaling pathways, is a 
major determinant of their resistance to apoptosis induc- 
tion by cytokines. Studies on pancreatic islet microorgan in 
toto will allow to experimentally confirm these results in a 
3D natural system. It also will permit to verify the hypoth- 
esis that increased TNFa synthesis by a cells during insu- 
litis protects them by activating survival pathways, while 
priming a deadly regulatory loop and causing |3 cells apop- 
tosis (see Additional file 13). 

Methods 

Cell culture and treatment with cytokines 

Mouse glucagonoma cell line aTCl-6 was obtained from the 
American Type Culture Collection (ATCC); it was cultured 
in complete Dulbeccos modified Eagles medium (DMEM, 
Sigma- Aldrich®, Saint Louis, MO, USA), as described [33]. 
Mouse insulinoma cell line pTCl also was from ATCC; cells 
were grown in DMEM with 25 mM glucose (Sigma- 
Aldrich®), supplemented with 2 mM L-Glutamine, 15% heat 
inactivated (HI) horse serum, 2.5% HI fetal bovine serum 
(FBS), 1% penicillin/streptomycin, in 95% humidified air-5% 
CO2 at 37° C. Cells were passaged once a week after trypsini- 
zation and replaced with new medium twice weekly. Treat- 
ment with cytokines (recombinant murine IL-1|3, specific 
activity 5 x 10^ U/mg, Preprotech, London, UK, UE; recom- 
binant murine IFN-y, specific activity 1 x 10^ U/mg, Prepro- 
tech; recombinant murine TNF-a, specific activity 1 x 10^ 
U/mg, Preprotech) was as described [34]. aTCl-6 cells (pas- 
sages 20-40) were seeded the day before treatment in 60 
mm dishes at a density of 4.5 x 10^ cells. 

Evaluation of apoptosis and necrosis 

Percentage of apoptotic or necrotic cells was assessed 
through fiow cytometry. Analysis was performed with a 
Beckman Coulter Epics XL-MCL fiow cytometer (Beckman 
Coulter®, Inc., Hialeah, FL, USA). Cells were collected, 
washed with phosphate-buffered saline (PBS), and stained 
with Annexin V-FITC/propidium iodide (PI) (Sigma- 
Aldrich®) in Annexin-V binding buffer, as specified by the 
manufacturer. EXP032 ADC Analysis^^ software package 
(Beckman Coulter®, Inc.) was used for data analysis. 

RNA extraction and HT quantitative RT-PCR 

Total RNA was extracted with Trizol (Lifetechnologies™, 
Foster-City, CA, USA), according to manufacturers 



instructions. RNA quantification was performed with a 
Qubit® Fluorometer (Lifetechnologies™). RNA for HT 
miRNA expression profiling was reverse transcribed into 
cDNA of 519 mouse-specific and 68 rat-specific miRNAs 
through Megaplex™ RT Rodent Primer Pool sets, and pre- 
amplified through Megaplex PreAmp Rodent Primer 
Pool sets (Lifetechnologies ). Resulting cDNAs were 
loaded onto TaqMan Low Density miRNA Arrays (TLDA) 
cards, according to manufacturers instructions. TLDA 
cards were run on ABI 7900HT Real Time PCR system. 
RNA for Single TaqMan miRNA Assays was reverse tran- 
scribed into miRNA-specific cDNA through TaqMan® 
MicroRNA Reverse Transcription Kit (Lifetechnologies ) 
and amplified using TaqMan® Universal Master Mix (Life- 
technologies™), according to manufacturers instruction. 
RNA for analysis of miRNA targets was reverse tran- 
scribed into cDNA through High Capacity RNA-to-cDNA 
Kit (Lifetechnologies™) and amplified through Fast SYBR® 
Green Master Mix (Lifetechnologies™), according to man- 
ufacturers instruction. Primer sequences are available 
upon request. 

Criteria for selecting downregulated or overexpressed 
miRNAs 

Data quality and quantification were computed using 
Real-Time Statminer® software (www.integromics.com) 
(Integromics, Granada, Spain). Multiple reference genes 
were used to normalize data: Genorm (integrated into Real- 
Time Statminer®) [35] and DataAssist™ (Lifetechnolo- 
gies™) softwares allowed to choose the best ones. Limma 
test (see below: Statistical analysis) was carried out by 
Real-Time Statminer® to assess statistically significant DE 
genes. DE miRNAs were ranked based on their p-values 
and adjusted p-values (Benjamini-Hochberg correction with 
False Discovery Rate, FDR, set at 5%). Relative quantities 
(RQ) of miRNAs between aTCl-6 and |3TC1, at steady 
state as between treated aTCl-6 cells and matched un- 
treated controls, were calculated according to 2'^^^^ 
method [36]. Values are reported as average fold changes of 
three independent biological replicates; RQ values < 1 were 
converted to negative fold changes by the formula: -1/RQ. 
The data files for each array are publicly available at the 
Gene Expression Omnibus (GEO) database repository 
(http://www.ncbi.nlm.nih.gov/geo/) (GSE42970). 

Transient transfection of aTC1-6 cells with mimics of miR- 
296-3p and miR-298-5p 

For transfection, aTCl-6 cells were plated into 24-well 
plates at a density of 4x10^ cells per well to obtain RNA, 
and into 100 mm dishes at a density of 1.15x10^ cells to 
obtain proteins. Transfections were performed using 
siPORT NeoFX (Lifetechnologies ) with 30 nM mimics 
of miR-296-3p/miR-298-5p/scrambled sequence (Pre- 
miR miRNA Precursor Molecules — Negative Control #1, 
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Lifetechnologies ). For each experiment, efficiency of 
transfection was measured through real-time PGR. 

In silico identification of miRNA targets 

VaUdated targets of DE miRNAs were retrieved from the 
Uterature and miRTarbase (release 2.5) [37]. Predictions 
were performed through starBase (release 2.1) (http:// 
starbase.sysu.edu.cn/). Among validated and predicted 
targets, only genes expressed in pancreatic cells (data 
from Beta Cell Gene Atlas, found at http://tldbase.org/ 
page/AtlasHome) and known to be functionally involved 
in cell survival or death were chosen for real-time PGR 
and western blot assays. Data on genomic position of 
genes encoding miRNAs and their assignment to specific 
families were from MiRBase (http://www.mirbase.org/). 

Western analysis 

Protein lysates and their quantification were obtained as 
previously described [38]. 50 (ig of total protein extract 
were loaded into 10% SDS polyacrylamide gel (Hoefer 
miniVE, GE Healthcare®, Amersham Place, Buckingham- 
shire, UK) and blotted to nitrocellulose membranes by iBlot 
Dry Blotting System (Lifetechnologies ). Membranes were 
probed with polyclonal antibodies to IGFlRp (Santa Gruz 
Biotechnology®, Inc., Dallas, TX, USA), p-IRS-1 (Santa 
Gruz Biotechnology®, Inc.), total-IRS-1 (EMD Millipore 
Gorporation®, Billerica, MA, USA), p-p44/p42 MAPK (Gell 
Signaling Technology®, Inc., Danvers, MA, USA), total 
p44/p42 MAPK (Gell Signaling Technology®, Inc.), TNFa 
(Gell Signaling Technology®, Inc.), using |3-actin (Sigma- 
Aldrich®) as loading control. Proteins were detected by 
using EGL™ Western Blotting Detection Reagents (GE 
Healthcare®). Densitometric analyses were performed by 
ImageJ software (http://rsbweb.nih.gov/ij/index.html) 

Prediction of transcription factors regulating expression 
of miRNAs 296-3p and 298-5p 

Matlnspector from Genomatix (http://www.genomatix. 
de/) was used to identify Transcription Factors Binding 
Sites (TFBS) and their corresponding Transcription Fac- 
tors (TFs) [39]. By using Multi Experiment Matrix (MEM) 
(http://biit.cs.ut.ee/mem) (selected collection: Afiymetrix 
GeneGhip Mouse Genome 430 2.0 [Mouse_430_2] plat- 
form, used to analyze 1546 datasets), TFs prediction was 
interpolated with data on statistically significant expres- 
sion correlation among TFs, which regulate DE miRNAs 
and their mRNA targets. Settings used in MEM are 
described in Additional file 14. 

Identification of CpG islands upstream the genes 
encoding Nespas, miR 296-3p and miR 298-5p 

GpG islands upstream genes encoding miR-296, miR-298, 
Nespas were identified through UGSG Genome Browser 
(http://genome.ucsc.edu/). 



Network analysis 

Biological networks comprising miRNAs, their predicted 
upstream regulators (TFs), their validated targets and 
first neighbours interactants, were generated by retriev- 
ing interactome data through MiMI Gytoscape plugin 
[40] and visualized by Gytoscape v. 2.8.1. Biological pro- 
cesses and pathways involving network nodes were ana- 
lyzed through the tool DAVID (http://david.abcc.ncifcrf. 
gov/) and BiNGO Gytoscape plugin [41] 

Statistical analysis 

P-values were calculated by applying different methods: 
Limma test [42], associated with Benjamini-Hochberg cor- 
rection for multiple comparison, was applied to identify DE 
miRNA genes between test and control samples in HT 
miRNA transcriptome analyses; Students ^-test, associated 
with Bonferroni correction method, was used to statistically 
analyze data from single TaqMan assays and to assess sig- 
nificantly different apoptotic levels between aTGl-6 treated 
with cytokines and their matched untreated controls; Tukey 
HSD post-hoc one-way ANOVA test was used to evaluate 
significant differences in apoptosis among different trans- 
fection experimental conditions. For analysis of correlation 
between the expression of miR-296-3p and miR-298-5p in 
aTGl-6, Pearson correlation coefficient was calculated. 
Finally, /^-square test was used to establish if miR-296-3p 
and miR-298-5p have more common targets than expected 
by chance; Fishers exact test was applied to evaluate the 
enrichment in specific gene ontologies. All statistical tests 
and correction methods, used to calculate p-values, are 
described throughout the text and figure legends. 

Nomenclature of genes and proteins 

Rules for official gene and protein symbols by the Inter- 
national Gommittee on Standardized Genetic Nomencla- 
ture for Mice were followed throughout [43] (http:// 
www.informatics.jax.org/mgihome/nomen/gene.shtml). 

Additional files 



Additional file 1: DE miRNAs between aTC1-6 and PTC1 at steady 
state. 

Additional file 2: Classification in families of miRNAs more 
abundantly or exclusively expressed in aTC1-6 respect to piCI at 
steady state. 

Additional file 3: Classification in families of miRNAs more 
abundantly or exclusively expressed in (JTCI with respect to aTC1-6 
at steady state. 

Additional file 4: aTC1-6 DE miRNAs after 24 and 48 h of treatment 
with cytokines. In bold are indicated significant DE miRNAs (adjusted p- 
values<0.05. Limma Test, Benjamini-Hochberg correction). MiRNAs are 
ranked based on increasing adjusted p-values (from the most to the least 
significant DE miRNAs). 

Additional file 5: Classification in families of DE miRNAs in aTC1-6 
after treatment with cytokines. 
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Additional file 6: Scatter plot showing correlation between miR- 
296-3p (x-axis) and miR-298-5p (y-axis) expression in aTC1-6, during 
a 6-12-24-48 h time-course experiment. For each time point DCt 
values of miR-296-3p and miR-298-5p were correlated, both from 
untreated and cytokines-treated aTCl-6 cells (r-value = 0.88, p-value = 
1.15e-08, Pearson's correlation test). Three independent biological 
replicates (n = 3) have been analyzed at each time point. 

Additional file 7: On scale representation of the genome segment 
comprising Nespas, miR-296, miR-298. CpG islands are indicated as 
red vertical lines; pre-miRNAs 296 and 298 are depicted as blue and 
green boxes, respectively; exons of noncoding RNA Nespas are shown as 
yellow boxes. Expression of a macro-noncoding RNA (precursor of miR- 
296, miR-298, Nespas) is predicted to be controlled by two groups of CpG 
islands (one comprising two CpG islands, from 17.5 to 18.8 kb upstream 
the first nucleotide of pre-miR-296; the other made of three CpG islands, 
from 30.1 to 33.6 Kb upstream the first nucleotide of pre-miR-296). 

Additional file 8: Validated and predicted targets of miR-296-3p 
and miR-298-5p. A selection of validated and predicted targets of miR- 
296-3p and miR-298-5p was chosen according to their involvement in 
apoptosis, cell cycle progression, cell differentiation and hormone 
secretion. 

Additional file 9: Interaction network among miRNAs more 
abundantly expressed at steady state in aTC1-6 with respect to 
PICI, their validated targets and first neighbours interactants. The 

network generated from validated targets and first neighbours 
interactants of 14 out of 50 miRNAs, more abundantly expressed at 
steady state in aTCl-6 with respect to pTCl, consisted of 560 nodes and 
8222 edges. MiRNAs are represented as fuchsia diamonds. This file can be 
opened and browsed through Cytoscape tool and its plugins (http:// 
www.cytoscape.org/). 

Additional file 10: Interaction network among miRNAs specifically 
expressed at steady state in aTC1-6 (a-miRNAs) with respect to 
PICI, their validated targets and first neighbours interactants. The 

network generated from validated targets and first neighbours 
interactants of 4 out of 23 a-miRNAs consisted of 1 1 7 nodes and 530 
edges. MiRNAs are represented as fuchsia diamonds. This file can be 
opened and browsed through Cytoscape tool and its plugins (http:// 
www.cytoscape.org/). 

Additional file 11: Interaction network among miRNAs more 
abundantly expressed at steady state in pTCI with respect to 
aTC1-6, their validated targets and first neighbours interactants. 

The network generated from validated targets and first neighbours 
interactants of 36 of 74 miRNAs more abundantly expressed at steady 
state in pTCl with respect to aTCl-6 consisted of 439 nodes and 2079 
edges. MiRNAs are represented as fuchsia diamonds. This file can be 
opened and browsed through Cytoscape tool and its plugins (http:// 
www.cytoscape.org/). 

Additional file 12: Interaction network among miRNAs specifically 
expressed at steady state in PTC1-6 (P-miRNAs) with respect to 
aTC1-6, their validated targets and first neighbours interactants. 

The network generated from validated targets and first neighbours 
interactants of 15 of 26 3-miRNAs consisted of 107 nodes and 132 edges. 
MiRNAs are represented as fuchsia diamonds. This file can be opened 
and browsed through Cytoscape tool and its plugins (http://www. 
cytoscape.org/). 

Additional file 13: Hypothetical model of regulation of miR-296-3p 
and miR-298-5p biomolecular activity in aTC1-6 at steady state 
(left) and after treatment with cytokines (right). 

Additional file 14: Additional Methods. 
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AT: After transfection; DE: Differentially expressed; DM: Diabetes mellitus; 
GLPl: Glucagon-like peptide 1; HT: High-throughput; IFN-y: Intereferon 
gamma; IL-1[3: Interleukin 1 beta; miRNA: microRNA; Ngn3: Neurogenin 3; 
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1 diabetes mellitus; T2DM: Type 2 diabetes mellitus; TLDA: TaqMan Low 
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